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The Rosenfeld functional provides excellent results for the prediction of the fluid phase of hard 
convex particle systems but fails beyond the freezing point. The reason for this limitation is the 
neglect of orientational and distance correlations beyond the particle diameter. In the current ar- 
ticle we resolve this restriction and generalize the fundamental measure theory to an expansion in 
intersection centers. It is shown that the intersection probability of particle systems is described by 
an algebra, represented by Rosenfeld's weight functions. For subdiagrams of intersection networks 
we derive vertex functions that provide the building blocks for the free energy functional. Their ap- 
plication is illustrated by deriving the Rosenfeld functional and its leading correction which is exact 
in the third virial order. Furthermore, the methods are used to derive an approximate functional 
for the infinite sum over Mayer ring diagrams. Comparing this result to the White Bear mark II 
functional, we find general agreement between both results. 
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I. INTRODUCTION 

Historically, density functional theory (DFT) for classi- 
cal particles has been investigated in parallel to its corre- 
sponding application in quantum mechanics [1] , although 
with less impact on our understanding of the underly- 
ing physics. One reason is a practical one, as Newton's 
equation is much more accessible by numerical methods 
than Schrodinger's. The expensive and difficult problem 
of constructing a suitable functional is therefore more 
profitable in the quantum case as for classical systems 
[2]. Nevertheless, the analytical form of a free energy 
functional contains rich information about the physical 
system that would otherwise be difficult to obtain from 
computer experiments alone. The construction of a clas- 
sical density functional is therefore of great interest from 
a theoretical point of view. 

An important step forward was the development of the 
fundamental measure theory (FMT) for hard convex par- 
ticles by Rosenfeld [MI], generalizing the semi-heuristic 
scaled-particle theory of Reiss, Frisch, and Lebowitz [5]. 
Starting from the observation that Mayer's / function is 
decomposable into a pairwise convolute of weight func- 
tions [SHU]) the free energy functional of uncorrelated 
particles is the sum of three contributions, constrained 
by the scaling behavior of the weight functions. This 
functional and its corrections by Rosenfeld and Tarazona 
[T2l [TB] proves to be of surprising accuracy. Compared 
to computer simulations, the phase diagrams for spheres, 
cylinders, rods, and their mixtures are in excellent agree- 
ment in the fluid region |14H19j and the direct vicinity of 
the freezing point [50]. On the other hand, the functional 
fails for higher particle densities. The reason for this 



* stephan.korden@rwth-aachen.de 



shortcoming is the missing correlation between orienta- 
tions and distances beyond the two-particle system. This 
restriction reflects the underlying Percus-Yevick approx- 
imation which breaks down for highly correlated particle 
configurations such as crystalline structures. 

Several important approaches have been made to ana- 
lyze and improve the Rosenfeld functional. A central step 
in this direction is the geometrically motivated correction 
term introduced by Tarazona |12j . Whereas a different 
approach compares simulation data to the structure of 
the functional, resulting in the White Bear functionals 
[211 [22]. A different strategy was followed by Leithall 
and Schmidt [23] by introducing a diagrammatic formu- 
lation relating the functional to a degeneration of Mayer 
diagrams. 

In a recent article, we started to investigate and clarify 
the mathematical origin of the local splitting of Mayer's 
/ function [2H US] . First it was shown that the kinematic 
formula of integral geometry, developed by Blaschke, 
Santalo, and Chern [26H30] corresponds to Rosenfeld's 
decomposition of the second virial integrand. More gen- 
erally, the intersection probability of any number of par- 
ticles, with a common intersection center, is determined 
by the Euler form and factorizes into a convolute of lo- 
cal densities. This result not only allowed the derivation 
of Rosenfeld's functional from first principles and with- 
out reference to the semi-heuristic scaled-particle theory, 
but also related the approach to Mayer's virial expansion 

In the current article, this connection between FMT 
and the virial expansion will be further investigated. It 
will be shown how to derive higher order terms of the 
free energy functional, extending the methods of [25J . As 
a first example, the leading order correction will be ob- 
tained, which resolves the angular degeneration of the 
Rosenfeld functional and clarifies how to correct the dis- 
tance correlation beyond one particle diameter and there- 



fore exceed the Percus-Yevick approximation. 

In order to keep the article self-contained, the algebraic 
results of [25] are repeated and refined in section [TT] and 
generalized to an intersection algebra in Sec. |III| Deriv- 
ing its representation in vertex functions in Sec. |IV| we 
calculate the leading correction to Rosenfcld's functional 
in Sec. [V] and finally compare the result to the White 
Bear mark II functional. 



II. REVIEW OF INTERSECTION NETWORKS 
AND THEIR EULER FORMS 

The recent investigation [25] has not only shown how to 
derive Rosenfeld's functional from the virial expansion, 
but also proposed a natural generalization to higher order 
corrections, where the methods developed for the leading 
order will also be of relevance for all further correction 
terms. We will therefore first give a short summary of the 
most relevant results obtained so far, including the graph- 
ical representation of intersection networks, the algebraic 
decoupling of the Euler form into weight functions, and 
their resummation into a generating function. 

The understanding of hard particle physics begins with 
the observation that the intersection probability of par- 
ticles, overlapping in at least one common point, around 
which the particles can freely rotate and translate, is de- 
termined by the Euler form. This central result of inte- 
gral geometry has been derived for two-particle domains 
by Blaschke, Santalo, and Chern [2TS1 - I5U] and further ex- 
tended to an arbitrary number of particles in [24 . The 
decoupling of the second virial integral into Rosenfeld's 
weight functions is therefore only a specific example of a 
more general relationship between differential geometry 
and the local Euclidean group ISO (3) of rotations and 
translations. 

With the exception of the second virial integral, the 
Euler form does not determine Mayer integrals exactly. 
However, completely connected Mayer clusters can be ap- 
proximated by the intersection probability of particles 
that intersect in at least one common point. The Euler 
form determines therefore an essential part of these im- 
portant Mayer integrals, but at the same calculational 
costs as the second virial integral itself. Nevertheless, 
completely contracted diagrams are only the leading or- 
der in an expansion of the free energy functional in 
the number of intersection centers, which constitute the 
smallest unit of an intersection network. Because of their 
importance, we introduced the name "stack" and "uni- 
versal stack" in [35J, defined by: 
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0< 



ust = st A 



(1) 



fc=2 



for the i = 1, . . . , k particle domains Di intersecting in at 
least one common point. 

The second virial integrand is equivalent to Mayer's / 
function. This guarantees an exact relationship between 
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FIG. 1. Mayer clusters and intersecting diagrams provide 
identical representations for pairwise intersecting particles, 
with nodes and edges interchanged. This is shown for the 
two diagrams in the Mayer (left), particle (middle), and in- 
tersection (right) representation. 
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FIG. 2. Integrals of intersection diagrams are approximated 
by consecutive contraction of their intersection centers. The 
example shows the three 4-particle Mayer clusters, their dual 
intersection diagrams, and their allowed contractions defined 
by the Mayer clusters. 



Mayer diagrams and their representation as intersection 
networks of only pairwise overlapping particles. For the 
simplest cluster diagrams these can be illustrated as 2- 
dimcnsional drawings. However, to simplify the graphi- 
cal representation, we also introduced "intersection dia- 
grams" in |25| . where particles are reduced to lines and 
intersection centers indicated by edge joints. An example 
with all three different types of representations is shown 
in Fig. [I] Approximations of these diagrams are derived 
by the successive contraction of intersection centers as 
shown in Fig. [2] for the four particle cluster diagrams. 

Intersection diagrams can be classified by their number 
of intersection centers h and internal loops g. Taking this 
into account, the excess free-energy functional density is 
the infinite sum 



=a.h= 



(2) 



where each element <& g ,h corresponds itself to an infinite 
set of diagrams. The leading element, ^o.ii is presented 
in Fig. [3j and provides the graphical representation of 
Rosenfeld's functional as the intersection probability of 
the universal stack. Observe that similar diagrams have 
been used in [53] to relate the virial expansion to FMT. 

The intersection probability for each stack St a, is de- 
termined by the Euler form A"(<9StA;), integrated over the 
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FIG. 3. The starfish like diagrams follow from the completely 
connected Mayer clusters by maximally contraction of their 
intersection centers. Resummation of all these diagrams is in- 
dicated by the crossed circle which corresponds to the Rosen- 
feld functional. 



intersection domain f a £ Stk and averaged over the po- 
sitions and rotations of each individual particle 



T(D) := {7 = (f,Q)\Fe D,Q£ SO(3)} 
dji := d 3 r l d 3 fli . 



(3) 



To keep the notation for the coordinates of particles and 
their intersection centers apart, their indices will be la- 
beled by the characters: 



a, b, c, . . . : indices of intersection domains 
i, j, k, . . . : indices of particle domains . 



(4) 



Suitably normalized [35] and multiplied by the single 
particle densities pi , the Euler form determines the virial 
coefficient at 0-loop order 
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if(c>St fe _i)<5(nr )d 3 r Q 

St k xr(D 2 x...xD k ) (5) 

x pi (71) . . . pk(jk) dj2 A . . . A d-fk , 



taking into account the symmetry factor Gk-i = 1 of 
the fully connected Mayer diagrams. The delta- function 
S(nf a ) of the scalar-product of the normal and position 
vector restricts the integrand to the surface dD and fol- 
lows the calculational rules introduced in Appendix B of 

The Euler form is a linear functional on the boundary 
of manifolds, which vanishes for any odd-dimensional do- 
main. The derivation of K(dStk) for 3-dimensional parti- 
cles thus reduces to only three contributions. Introducing 
the surface £ of a domain D and the notations: 



dD , s" = En...nE. 



(6) 



n— times 



the boundary of a stack of k identical particles reduces 
to 

dSt k = k £ n St fc _i + k(k - 1) £ 2 n St fe _ 2 

+ /s(/ J -i)(fc-2)i] 3 nst, £ _3 

as the intersection probability between the set of points 
£ 3 and a further surface element K(T, 3 n £) = is zero. 
The only non- vanishing terms of the Euler form are there- 
fore: 



X(£) = u x 

K{Y?) = C a ^LO ai uj a2 

K (£ 3 ) - C^'^u^u^u^ 



(8) 



with the C-matrices C" 1 " 2 , C""" 2 " 3 explicitly derived 
in [25]. As first shown by Chern [28], these tensors are 
independent of the particle geometry and solely defined 
by the Euler form and the dimensions of the particles and 
their embedding space. 

To this list of algebraic relations it is useful to add a 
further one, which follows from the tensorial density of 
the integrand ([5]): 



K(T, n n Stfc) = K(Yi a ) {u) v ) 1 



(9) 



For 3-dimensional convex particles it has been shown 
in [35J that the infinite dimensional basis set of weight 
functions can be grouped into five classes: 
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= Q(f a - fi) , 



(10) 



corresponding to the Euler-characteristic x of the Gauss 
curvature kq, the mean curvature R, the curvature dif- 
ference or tangential curvature A, the surface a, and the 
particle volume v. Each of these geometric terms is taken 
at the intersection center r a with respect to their absolute 
position in the embedding space fi. As a consequence of 
the non-algebraic splitting of the scalar Euler form, the 
weight functions also depend on the L-fold tensor prod- 
uct of the normal vector n for L £ Nq. 

The weight function ui v of the particle volume v plays 
a special role. It is not part of the curvature dependent 
Euler form Q but constrains the integration domain in 
^ from the embedding space M 3 to the particle volume. 
The relation ([9| is therefore a formal one. Nonetheless, 
it is useful to include it to the set of algebraic relations 
^ and to introduce two different indices for the weight 
functions 



A, B, C . . . £ {v, x, kL, <jL} 
a,fi,j...£ {x,KL,aL} , 



(11) 



indicating if v is included or not. The characters are cho- 
sen such that to each intersection center f a G Stfc the 
family of indices (a, A\, ax, A 2 , a>2, ■ ■ •), (b, Bi, j3\, . . .), 
(c, Ci, 71, . . .), ...is assigned, providing an intuitive re- 
lation between weight indices and intersection points. 

It is a special property of the 0-loop order ^ and its 
single intersection center that each particle domain Di 
is related to a single weight function uj\. Only in this 
special case, it is possible to combine them into the 1- 
point density 
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(12) 
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introduced by Rosenfeld [3 . Whereas higher loop orders 
require the definition of the "fc-point function" 



(13) 



for k disjunct intersection centers r a at particle domain 
fi G Di and its corresponding "fc-point density" 



k ( r ai i 



M k 



i) drfi , 



(14) 



as introduced in [55], generalizing Eq. (12 1 and 
Wertheim's 2-point measure [3"2T - f3"5] . 

Due to the coupling of the particle density to the 
weight functions, the free energy F is no longer a func- 
tional of pi alone. Instead, F now depends on the new 
variable n v . As has been shown in [35], the weight 
function u> v acts as the neutral element under remov- 
ing Stfe — > Stfc_i or adding St k — > Stfc+i a particle and 
thus shifting the Euler form ^ by one factor of u> v . The 
corresponding shifts for the boundary of the stack ^ or 
its Euler form 

K(dSt k ) = k K{T) u k v - x + fc(fc - 1) K{Y?) w k v - 2 

+ k{k-l){k-2)K(Y 3 )u k - 3 ' ''' 

are generated by integration and differentiation with re- 
spect to to v : 



J K(dSt k )dw v = ^L^K(dSt k+1 ) 



SK(dSt k ) 
5oj v 



(16) 



kK(dSt k -! 



and thus relate intersection integrals |5| for more than 
three particles k > 3. These operations allow to translate 
the virial expansion in the particle density representation 
p to that in the weight density n v . To see this, consider 
the virial expansion of the chemical potential [31 , written 
for constant p and /3 _1 = keT: 



Pk-i = 
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fip = (3p id + y]fik-il 
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(17) 
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depending on the embedding volume V and the symme- 
try coefficient a k . 

From this derives the free energy potential by adding 
one further particle, realized as the integral over p: 



/3F = /3F id + (3F e: 



oo „ 

f> F - = E / u Pk - ipk 

fc=2 J 



(18) 



In the representation of weight densities, the last step 
corresponds to the shift Stfc_i — > St k in the particle stack, 
which for k > 3 particles is realized as the integral over 



The operations ( 16 1 therefore apply to the free energy 



functional, defined as the generating function of cluster 
integrals. With the functional derivative 



Sn A (r a ) 
Sn B (f b ) 



= 5ab S{r a - n) 



(19) 



the free energy F is related to the free energy density $ 
by the integral 



0F = 



4(,5F W„)= / 'a^A. 



Sn v (r a ) 
®(r a ) d 3 r a 



Sn v (f a ) 
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§ g ,h{r ai ,...r ah )d 3 r ai ...d 3 r ah 



(20) 



which reduces to the expansion in intersection diagrams 
by Eq. Comparing this result to the virial represen- 
tation of the free energy ( 18 ) yields a relation between the 



density functionals <& g .h and their corresponding cluster 
densities 



f °° 1 

i>([rivt],r a ) = / ^ r l3 k -i([n A ],f a )dn v 

J k=2 



(21) 
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that is uniquely defined up to an integration constant. 
Its value has been determined in |25| and corresponds to 
the formal definition of a single particle virial coefficient 



A> = ft 



(04) 



(22) 



The final step in proving the equivalence of the inter- 
section probability of the universal stack and Rosenfeld's 
functional consists in deriving the first element $o,i of 
( 20 1 . Inserting the algebraic representations p| , ^ into 
the Euler form ( 15 1 and rewriting the cluster integral ^ 
in 1-point densities 



, M I, ,k— 1^2— ik 

W Y (^V ) 



+ (k - 1) r ia! «(^ 2 ) i3 -" 

+(k - l)(fc - 2) C QlQ2Q3 « 2 ^ 3 (^- 3 ) 44 - 1 ' 
x p h ... p lk dj i2 . . . d-y ik d 3 r a 



(23) 



n x n k ^ + (k - 1) C aia2 n ai n a2 nl 



k-2 



+ (k - l)(fc - 2) C aia2a3 n ai n a2 n a3 n 



k-3 



d 3 r n 



5 



yields the decoupled integral for a stack of order k. 
Adding up all cluster integrals and integrating over the 
packing density n v 

*0,i([»a],»*«) = / E V^([nAlr a )dn v (r a ) (24) 
^ 7 1 



n 

1 — n 
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(l-n v ) 2 



-2C* 
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reproduces Rosenfeld's functional 
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and identifies it as the 0-loop order of the expansion in 
intersection centers. 

The simple structure of this result is explained by the 
single intersection point, as shown in Fig. [3] However, 
higher order intersection diagrams, as shown in Fig. [2] 
not only incorporate further intersection points but also 
loop constraints that create corrections to the direct cor- 
relation function at distances larger than the particle di- 
ameter. Going beyond the Percus-Yevick approximation 
requires therefore the introduction of additional mathe- 
matical tools. 



III. THE INTERSECTION ALGEBRA 

Intersection diagrams were introduced in |25j to visu- 
alize the approximation scheme of FMT and to associate 
the virial expansion of the free energy to the Euler form 
and thus to Rosenfeld's functional. They also gave a 
first intuitive understanding of higher order corrections as 
partially contracted diagrams. Figures Fig. [I] and Fig. [2] 
demonstrate how they can be obtained by graphical con- 
struction. For more complex diagrams, however, this ap- 
proach becomes unwieldy and algebraic rules for their 
construction and manipulation are more useful. 

As an example, we will first reconsider the intersection 
diagrams of Fig. [I] and Fig. [2] in detail in the first para- 
graph |III B| and then generalize the results to diagrams 



of arbitrary degree of contraction in III C 



A. Intersection diagrams 

The Euler form provides a unique identity between 
Mayer's / function and the weight functions 



fij(n-fj) 

-S 



C A ^u l Al (f a - fi)J A2 (r a - r 3 ) d 3 r a (26) 



and relates the representation in particle r^, fj and inter- 
section coordinates r a . Thus any sequence of / functions 
/12/13/14 ■ • ■ , multiplied by their particle densities, can 
uniquely be rewritten in fc-point densities ( 14 ) as a func- 
tion of their intersection coordinates. 

The simplest diagram, apart from the second virial 
cluster, is the triangle graph, which has been discussed 
by Wertheim 32-35] and in [25]. Using the indices of 
Fig. [4)3) , its corresponding integral in weight functions 

1 



(27) 



x C B ^uj 2 B2 {r b ~ r 2 )uj 3 B3 (r b - r 3 ) 
x d 3 r a d 3 r b d 3 r c p 1 (71)^2 (72)^3 (73) djid^d^ 



is a functional of the particle densities. The same inte- 
gral in 2-point densities (14 1 allows the more compact 
notation 
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_qA 1 A 2 (jB 2 B 3 



c c 3 c 



(28) 



x J n Al cJr ac )n A2 B 2 (r ab )n B3 c 3 (n c )d 3 r ab d 3 r bc d 3 r ca 

using the distance vectors f ab = f b — f a . 

In the following discussion neither the dependence on 
the particle densities nor the loop constraints will be of 
relevance, so that the identity ( 26 1 can be written in the 
simplified form 



fi2(A)=C AlA2 uj Al uj A2 , 

using the index combination 

Ai : weight index A at intersection 
domain r a of particle i 



(29) 



(30) 



and omitting any reference to the integration over the in- 
tersection coordinates. The index A has now two mean- 
ings. On the one hand it indicates the intersection center 
A, on the other hand it also numbers the link of the Mayer 
cluster fij{A), as seen in Fig. [4p 



The two representations (|27j) and ( 28 ) of the three par 



ticle Mayer cluster of Fig. Wp) can now be written in the 
more convenient form: 

/i 2 (A)/ 23 (B)/3i(C) (31) 

=C mm u a ^a 2 ■ C B2B3 u B2 ojb 3 ■ C C3Ci uj C3 lu Ci 

= C AlA<2 C B ^ C c ^ { UAl uj Ci ) {u A2 u B2 ) (u Ba u Ca ) 

with the definition of the 2-point function ( p~3| ) indicated 
by parenthesis. Correspondingly, example Fig. Hp) is the 
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FIG. 4. Examples for the index conventions for Mayer and 
intersection diagrams. For pairwise intersecting particles, the 
two representations are dual to each other with the role of 
nodes and edges interchanged. Upper case letters indicate 
Mayer bonds or intersection centers. Whereas particle num- 
bers represent root points or particle lines. 
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FIG. 5. Substitution rules for Mayer diagrams: a) The Mayer 
vertex of particle i with k attached bonds corresponds to a 
fc-point function and b) the Mayer bond linking particles i 
and j maps to a C-matrix. 



product of 3-point and 2-point functions 
fi2(A)f 23 (E)f 34 (C)f 41 (D)f 24 (B) 

^ < jA 1 A 2( jE 2 E 3 (jC 3 C i( jD 1 D i< jB 2 B i (0,2) 

x (uj Ai uj Di )(w A2 ujb 2 ^e 2 )(^c 3 ^e 3 )(^c^b 4 0Jd 4 ) ■ 



Comparing the three examples (29), ( [31] ), and (32) 
to their corresponding Mayer diagrams in Fig. [4] reveals 
a simple building rule for virial integrals represented in 
weight functions. Define the "Mayer vertex" as the node 
of a Mayer diagram with its attached edges, as shown 
in Fig. [5^,). The cluster integral can then be read off 
directly from the Mayer graph by the substitution 



Mayer vertex — > fc-point function 
Mayer bond — > C-matrix 



(33) 



multiplied by the particle densities, loop constraints, 
integral measures, and the symmetry factor. For the 
generic case of pairwise intersecting particles, this sim- 
ple replacement defines a unique mapping between the 
representation of Mayer and intersection diagrams. 

Graphically, the building rule ( 33 1 corresponds to an 
exchange of nodes and edges, as shown in Fig. |4j which 
in terms of graph theory defines the " dual graph " [36 of 
the Mayer diagram. However, this bijection still provides 
no simplification of the integral, and the evaluation of the 
dual diagram is as complicated and unmanageable as it 
is for the original virial cluster. The next step considers 
therefore the systematic approximation of intersection di- 
agrams of pairwise intersecting particles. 




b) D 



FIG. 6. The approximation scheme of FMT is based on the 
successive contraction of intersection points. Reducible dia- 
grams are intermediate steps with two particles intersecting 
more than once. Further contractions finally result in irre- 
ducible graphs. 



B. Contraction rules 

It has already been shown in [25] that the structure of 
intersection diagrams can be simplified by moving some 
of their intersection centers into a single one. This pro- 
cess of "contraction" increases the rotational and trans- 
lational degrees of freedom over which the statistical sys- 
tem is averaged. It therefore reduces the complexity of 
the virial integrals but also coarsen their spatial resolu- 
tion. For the most simplest cases, as the clusters of Fig. [2] 
and Fig. [3] the contraction of diagrams can be done by 
hand. 

To derive an algebraic set of contraction relations, con- 
sider the pairwise contraction of the triangle diagram 
shown in Fig. [6^,). The three intersection centers are 
combined in two steps: First, C is shifted into B, leaving 
the point A invariant; then B is shifted into A. For each 
of the three diagrams of Fig. [6^,) we can now write down 
their corresponding virial integrand, using the splitting 
rules ([8| and ^ of the Euler form: 



B^A 



> C 



MM 



W Al U}A: 



■ C 



B2B3B1 



> c 



A X A 2 A 3 



u Al u; A2 u} A3 



lu B2 ujb 3 ^b 1 (34) 
(35) 



The first rule of the contraction operation can now 
be read off from (34) for the fusion of two intersection 
centers 

lim C B2B3 C C3Cl oj B2 ub 3 uc 3 u Ci 
C ^ B (36) 

= C B2BsBl U B2 UJB 3 U Bl ■ 

More clearly, the effect of the limit C — > B on the 
weight functions and C-matrices can be grouped into two 
classes. If the objects belong to different particles 



lim cjo„ = lub 
c^b J 

lim C C3Cl = C 

C^B 



3 ! 

B 3 B 



lim ujb 2 uj C3 = ub 2 Ub 3 
o — y -D 



(37) 



the indices of the intersection centers are simply renamed. 
The same operation for identical particle indices 



lim C B ^C C ^ 



C^B 



wc 3 



C 



B2B3B1 



UJ Ba 



(38) 
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yields a tensorial contraction of the C-matrices, aligned 
with the removal of one weight function. This operation 
is only defined in the combination of weight functions 
and C-matrices and cannot be split in the way of (37 1. 
Applying these rules to Eq. ( 34 ) 



lim C AlA2 C B2B3Bl lj Ai lu A2 lj B2 oj B3 lj Bi 
B ^ A (39) 

= C A lA2 A 3 u u u 



reproduces the graphically obtained result of ( 35 ) , pro- 



vided Eq. ( 38 1 is extended to two identical particle in- 



dices. This observation is readily generalized to the pair- 
wise contraction of n — k coincident particle indices 



lim c Al '" AkAk+1 '" A " Q B k+i—B n B n+1 . 



x u Ak 



C 



A-L...A k 



.A n ...A m . , 



■ u >B n 
■ ■ ■ U A„ 



(40) 



Thus the successive contraction of intersection centers 
generates higher rank C-matrices. However, from the 
splitting rules of the Euler form ^ and Q we know 
that the maximal rank of a C-matrix for 3-dimensional 
particles is at most 3 and that all further indices neces- 
sarily reduce to the index v of the particle volume. It is 
therefore natural to combine the two equations ^ 
into one 



K(dSt k ) = Ct A *~ A Ju Al ...u At 
and to define the generalized C-matrix 

c (A H A i2 Ai 3 A i4 ...A ik ) ._ c ( Ail A i2 A i3 



4: J 



(41) 



(42) 



where the parenthesis indicate the symmetrization of all 
particle indices. 



element u) v for three identical particles 



As an example, let us expand Eq. (35) in the neutral 

(43) 



+ C (a *i 



s 3 

.Wo 



c< 



3 uj x ujI + 6 C ai ° 2 uj ai uj a2 uj v + 6C Ql 



The result correctly reproduces the Euler form ( 15 ) for 
k = 3. 

The successive application of pairwise contractions on 
dual Mayer clusters generates a vast number of diagrams. 
However, some of them correspond to networks with mul- 
tiple intersections between particles, as shown in Fig. [6^,). 
The first step C — > B generates an intermediate diagram 
with the particles 1 and 2 intersecting twice in the cen- 
ters of A and B. But as one intersection point already 
determines the position and orientation of their parti- 
cles uniquely, this diagram is no allowed configuration. 
Whereas the next contraction, B — > A, resolves this am- 
biguity. 



Intermediate diagrams are identified as products of C- 
matrices with more than one common particle index. To 
distinguish these cases from admissible intersection dia- 
grams, we introduce the notation: 

( jA 1 ...A k A k + 1 ...A n QB k+1 ...B n B n+1 ...B m 

n — k = 1 : irreducible intersection (44) 
n — k > 1 : reducible intersection. 

Correspondingly, diagrams without reducible intersec- 
tions are referred to as "irreducible diagrams" and "re- 
ducible diagrams" otherwise. It follows from their defi- 
nition that any reducible intersection can be reduced to 
an irreducible one by further contractions. 

Another example is the Mayer diagram of Fig.|4j:). Its 
contractions can be either derived by ( 40 ) or read off from 
Fig.fi>) 



(J Al A 2 (jB 2 Bi fjDi Di qC 3 d (jE 2 E 3 

X U Al UA 2 U B2 U B4 UJ Dl UJ D4 U C3 U Ci UJ E2 UJE 3 
g ^C > ( jA 1 A 2( jB 2 B 4 . ( jD 1 D i( jC 2 C 3 C i 



X WA 1 WA 2 WB 2 WB 4 Wfl 1 WD 4 Wc 2 WC 3 WC 4 



A 1 A 2( jD 1 D lir ,B 2 B 3 B, 



(45) 



(46) 



X Wyt; W A 2 "Di wd 4 WB 2 Wfi 3 WB 4 - 

The first contraction, E — > C, yields again an interme- 



diate diagram ( 45 ) , reducible in the particle numbers 2 



and 4, which is then transformed into an irreducible one 
by shifting C — > B. The result is the highest possible ap- 
proximation of the Mayer cluster of Fig. |lj;) . No further 
contraction is possible as the particles 1 and 3 do not 
interact directly. This can be seen either from the Mayer 
diagram, where the corresponding /13 is missing, or di- 
rectly from the intersection diagram of Fig. [6)3) . There- 
fore, apart from irreducibility, the intersection diagrams 
obtained by pairwise contractions also have to be com- 
patible to the bonding relations of its Mayer graph. 

Again, compatibility of a diagram can be directly read 
off from its corresponding Mayer cluster, as each of 
the generalized C-matrices belongs to a completely con- 
nected Mayer subdiagram: 

(j{Ai— A k ) ~ completely connected Mayer 



subdiagram of the particles 1, 



,k. 



(47) 



The diagram Fig. |4j;) from the previous example can 
therefore be contracted either to C B2BiB * or C BlBlBi , 
corresponding to their subdiagrams of particle indices 
(2,3,4) and (1,2,4). 

In the next subsection it will be shown that any dia- 
gram can be decomposed into its maximally connected 
subgraphs. In most cases this splitting is uniquely de- 
fined. But the current example is one of the exceptional 
cases, which can be split in at least two different ways, 
corresponding to a global Z 2 symmetry. For such dia- 
grams, it is necessary to count their multiplicity of con- 
tractions, indicated by 



contraction multiplicity. 



(48) 



For the graph of Fig. Et) , the multiplicity is therefore 
m(l,0,0) = 2 and m(ni, n,2, = 1 else. (49) 

C. Contraction rules and their algebra 



The Euler form induces a real, linear operation on the 
set of Mayer clusters 

K(xiT ni + x 2 T n2 ) = xiK(T ni ) + x 2 K(T n2 ) 

K(T n ) = Y[K(T^) (55) 



The contraction rules are local mappings on the set of 
Mayer and intersection diagrams. In order to describe 
their operations on general graphs, it is practical to in- 
troduce a suitable notation for both types of representa- 
tions. 

As only completely connected subdiagrams can be con- 
tracted into single intersection centers, they take up the 
position of prime elements in the set of Mayer clusters 
and intersection diagrams. Let us therefore introduce 
the notation: 

: completely connected Mayer subdiagram of n 
particles and external bonds grouped in the 
partition A. 

Due to the permutation symmetry of the particle indices 
of completely connected diagrams, it is sufficient to group 
the external links into a partition table A. For example, 
A = [(ABC), (DEFG)} assigns the external links ABC 
to subdiagram 1 and DEFG to subdiagram 2. 

Any Mayer cluster can now be represented as a product 
of prime subdiagrams. For example, Fig. |4j;) allows the 
two decompositions: 



(ABD)CE r CE _ j,(BCE)AD r AD 



r? 



(50) 



which directly translates to the splitting of the Euler form 



(ABD)CE r CE\ _ 



{ABD)CE-s 



)K{T C 1 E ) (51) 



This notation is far more compact than the representa- 
tion in weight functions (45), (46). 



Analogously, intersection diagrams can be split at each 
intersection center along their particle lines: 



f A 
r 1 



intersection subdiagram of center A and 
partition A of external particle lines. 



(52) 



As an example, the intersection diagram of Fig. |4p) fac- 
torizes into the prime elements 



r 12 t^23 t^24t^34 t^41 
L A L E L B L C L D 



(53) 



In summary, the derivation and approximation of in- 
tersection diagrams reduces to a simple set of operations 
and constraints. In combination with the splitting rules 
of the Euler form and its linearity, we now define the 
"intersection algebra": 

Definition III.l Let K be the Euler form and T n ,b G T 
an element of the set of Mayer star-clusters T with n 
particles, factorizing into the prime subdiagrams 



E ; 



(54) 



for X\,x 2 G K, defining the "intersection algebra " (K, T). 
The splitting relation 



<^A 1 ^>A 2 ■ ■ ■ UA n 



(56) 



induces a representation on T in weight functions. Inter- 
section centers can be combined by pairwise contractions 

C:A^A', |A| > |A'| , (57) 

reducing the length of the partition |A| by at least one. 

The intersection algebra changes the focus from differ- 
ential geometry to the representation theory of the sym- 
metric group |37[ 138) , with the Euler form ( 56 ) relating 



the partition table A of an intersection diagram to the 
ring of symmetric polynomials. To prove that the con- 
traction operation respects this representation, apply the 
relation 

00 00 00 n 

(»•'•'-»•'•'' = E(E a ^-»)^ ( 58 ) 

i=0 j=0 n=0 k=0 

to the product of two generating functions, contracted in 
the first particle index i\\ 



lim (V C {M 

B^A ^ 



A 2 ...A n ),,h 



J A^A 2 



B\ B; 



UA n Pi 2 ■■■Pi n ) 

Pjh) Pil 
"LPh - --Pin . (59) 



which again corresponds to the Euler form of a new in- 
tersection diagram. 



The ring structure of the polynomial ( 56 ) greatly sim 



plifies the following derivation of the symmetry factors 
as well as the construction of vertex functions. 



IV. INTERSECTION VERTICES AND VERTEX 
FUNCTIONS 



The mapping ( 55 ) uniquely defines the splitting of any 



intersection diagram into its Euler forms. In principal, 
this is all one needs to derive higher order corrections of 
the free energy. However, much of this approach can be 
simplified by the resummation of subdiagrams. In the 
following two subsections, the free-energy representation 
in intersection centers will be generalized. First, it will 
show in |IV A| that the functional splits locally at each 
intersection center into vertex functions, whose analytical 
form will be derived in IIVBI 



A. The local splitting of the FMT functionals 

In the notation of the completely connected intersec- 
tion diagrams (52), the Rosenfcld functional is repre- 



sentable as the infinite sum 

r 00 1 

3 ,i(fA) = /E-^n 



• Pi n dn v (60) 



over the subclass of diagrams with only one intersection 
center. An even more concise notation can be obtained 



using the linearity of the Euler form ( 55 ) and considering 
the weighted sum over diagrams 



00 

E 

71=1 



1 



(61) 



symmetrizised over the external particle indices. This 
shortened notation provides a convenient representation 
for the discussion of diagrammatic resummation. 

But the Euler form is only one aspect in the deriva- 
tion of the free energy. In the following we will also need 
to generalize the expansion of the functional in intersec- 
tion centers ( 20 ) and to determine their combinatorial 



prefactors (21 ) of the virial contributions. 

Let us first focus on the expansion of the functional 
itself. Explicitly written out up to three intersection cen- 
ters 



0F = 



H/3F) 
5n v (f a 
<5 3 (/3F) 



Sn v (f a ) 



9=0 



8n v (f a )5n v (f b )Sn v (f c ) 



5n v (f a )5n v (f b )dn v (f c ) 



= $ 



$1,3 d 3 r a d 3 r b d 3 r c 



+ . . 



(62) 



it reproduces the integral representation of the Rosen- 
feld functional ( 24 ) at first order and yields the next to 
leading order 



( 00 1 

« I — 1 



(63) 

x dn v (f a )dn v {f b )dn v (r c ) , 



fol- 



which parallels the structure of (|27j for k = 3. 

From the local property of the Euler form (55) 
lows that the virial contribution factorizes into a prod- 
uct of three polynomials, each depending on one of 
n v (r a ),n v (r b ),n v (r c ). Consequently, the integration of 



( 63 1 factorizes likewise and can be executed for each in- 
tersection center individually. The same argument ap- 



plies of course to all further terms of the expansion ( 62 ) . 
The free energy functional decouples into a product of lo- 
cal functionals for each intersection center, coupled only 
by the loop constraints and the numerical prefactors of 
the virial integrals. 






v ; 



V 



FIG. 7. Two examples of the class of triangular diagrams, 
represented as Mayer graph, maximally contracted intersec- 
tion diagram, and as weighted graph. Their automorphism 
groups are: a) f(4,3,2): Ef x S 4 x S3 x S 2 and b) f(3,3,3): 

Gl2- 



Given the irregular structure of the symmetry factors 
cr(r) entering ( |17[ ) for different diagrams, it is nontrivial 
that such a splitting of the functionals § g ,h should ex- 
ists. A general proof of this hypothesis would require a 
classification of the automorphism groups of Mayer clus- 
ters under relabeling, which to the best of our knowledge 
is unknown. We will therefore focus on the leading cor- 
rection term of the Rosenfeld functional and explicitly 
derive ( 63 ) in the following sections. 

The diagrams entering $ 13 have three intersection cen- 
ters grouped into a ring and thus have a triangular sub- 
structure. Two examples are shown in Fig. [7j using the 
three representations as a Mayer graph, as the maximally 
contracted intersection diagram, and as a weighted dia- 
gram with the number of external lines as index. Using 
the notation of intersection subdiagrams (52), any such 
element can be written as 



r(m, n 2 ,n 3 ) := T 



lllz]l Jni p*2*3fcl— kn 2 p*3*lil— In 

B L C 



(r^ 1 ) iii2 (r^ 2 ) i2i3 (Tq 3 ) i3ii 



(64) 



with the paired indices ii, 12, £3 defining the backbone of 
the graph and j, k, I indicating its external lines. 

As shown in appendix [XJ the symmetry factor of a 
Mayer diagram T n of n particles is determined by the 
quotient 



<r(r n ) 



|Aut(r„)| 



(65) 



of the dimensions of the symmetric group S n and the 
automorphism group Aut(r n ) of the graph [39] - This re- 
mains true for the dual diagram T as the invariance group 
is independent of the representation of T and therefore 
unaffected by contractions 



Aut(r„) = Aut(r„) 



(66) 



However, our definition of the particle stack ([!]) and its 
representation as symmetrized C-matrix ( |42[ ) already in- 
cludes the invariance group of the completely connected 



10 



subdiagrams. It is therefore necessary to define an "ef- 
fective" symmetry factor for intersection diagrams a{T n ) 
without the invariance group of external particle-lines. 

As an example, consider the graph T(3,3, 3), shown 
in Fig. [Tja) . Each of the three intersection centers has 5 
particle-lines attached, of which 2 are fixed by the tri- 
angular backbone diagram, whereas the remaining 3 ex- 
ternal lines are invariant under the permutation S3. In 
order to compensate for the symmetry of the C-matrices, 
the group S3 x S3 x S3 = (S3) 3 has to be factored out 
from Aut(f(3,3,3)) = G 12 : 



G 12 /(S 3 ) 3 = D 6 



(67) 



The quotient group is therefore the automorphism group 
of the weighted graph shown in Fig. (7|d). As demon- 
strated in appendix [AJ this result applies to any dia- 
gram T(n,n,n), independent of the number of external 
particle-lines. 

Generally, it is far easier to derive the reduces invari- 
ance group from the weighted diagrams, where the ex- 
ternal particle-lines of the intersection graph have been 
replaced by their number as weight index. Using the 
results of Tab. |TTJ the "reduced automorphism groups" 
of triangular diagrams r(m, ^2,^3) consists of only four 
cases: 



m=n 2 = n 3 
n>l = n 2 = 0, 71,3 = 1 



El 



E x x Z 2 
S 2 x Z 2 . 



(68) 



Applied to the triangular diagrams (64), the partition 



table A = (Xa, As, Ac) consists of the three elements 



A.4 = ((11,12), (ji, ■ ■ • , Jm-2)) , 
As = {{i 2 ,h), (h, . . . , fc„ 2 _ 2 )) i 
Ac = ((h,h), (h, ■ ■ -,ln 3 -2)) , 



(71) 



from which follows the corresponding orthonormal set of 
unpaired indices 



A 



A = (il,--- J Jn 1 -2) , A^ = (fa, 

Ac = ■■■■> ^3-2) ■ 



, k n2 - 2 ) , 



(72) 



Using this notation, we define the effective symmetry 
factor of intersection diagrams: 



5(r„) :=^(r n )ni A ^l ! 

i 



(73) 



|Aut(r„) 



which determines the inverse of the dimension of the re- 
duced invariance group and therefore depends only on 
their equivalence classes of diagrams. Its values for the 
triangular diagrams are listed in Tab. [TTJ. 

This definition of the symmetry factor not only com- 
pensates for the symmetry of the C-matrix, but also in- 
cludes the normalization factor l/(n — 1)! of the virial 
coefficient (17) and the 1/n of its integral (18). All nu- 
merical prefactors of the free-energy virial expansion are 
thus included in a. 



with the exceptional diagram r(l,0, 0) first discussed in 
the context of the contraction rules in Sec. MI Bl 

The reduced invariance group identifies four equiva- 
lence classes of diagrams, independent of their particle 
numbers. This is an important property, as the final 
goal of determining the free energy functional $13 re- 
quires the resummation of all such diagrams. And here 
we see that this problem reduces to at most four different 
classes. 



Instead of the symmetry factor for virial integrals ( 65 1 , 



a depends only on the equivalence classes of diagrams 
and is independent of their particle numbers. In the rep- 
resentation of intersection diagrams (52), the number of 



external, unpaired particle-lines can be determined from 
the partition table A: 



Ai ' 



A = {X 1) ...X k } 



(69) 



Let (Xi)j € A, denote an individual element of the parti- 
tion Ai G A. The set of external lines is then character- 
ized as: 

A± := { (A,), I (Xi)j e A, , (A,-), g A fe , V i ? k} . (70) 



B. Vertex functions 

The new symmetry factor a does not affect the previ- 
ous derivation of the Rosenfeld functional ((23|, (25). Its 



numerical value for the starfish diagrams with unpaired 
particle-lines Ai = A^ 



a(VV 



1 



1 



(74) 



coincides with the corresponding symmetry factor of the 
Mayer diagram er(r^) = 1- The structure of the re- 



summation of 0-loop diagrams (61 1 remains therefore the 
same: 



QVV fi 1 ...i n = ff W p(H...g 



n=l 

00 



n ^ 



(75) 



(»!...*„) 

A 



n=l 



To get a first impression, how this result generalizes 
to diagrams with more than one intersection center, con- 
sider Fig. [8] The first element of the series of triangular 
diagrams is the exact third virial integral, which later 
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$1,3 = 



A+A + A + A_+ 



i 3 1 \ ii 



FIG. 8. Resummation of the triangular diagrams provides the 
first order correction to the Rosenfeld functional. With the 
intersection centers indicated by crossed circles, the functional 
is the product of three 2-vertex functions summed over the 
inner particle indices 



will replace the approximate term in the Rosenfeld func- 
tional. The next element is the contracted diagram of 
f (1,0,0), shown in Fig. §d) and obtained by attaching 
one additional external particle line to one of the three 
intersection centers A, B, C. Iterating this operation on 
the backbone diagram, it generates the series of $1,3, 
which can be resummed in the same way as the 0-loop 
diagrams, shown in Fig. [3J For each intersection center 
we thus obtain a "vertex function" A — > Vi r i 2 (A) depend- 
ing on the particle- indices ii, ta? *3- 

From an algebraic point of view, this factorization of 
the functional is a consequence of the decomposition of 
the triangular diagram (64 1 into subdiagrams. However, 



although uniquely defined for a given set of particle in- 
dices, it is invariant under permutation of the intersec- 
tion centers A,B,C. Instead of the single Mayer graph 
r(l,0,0), there exists three identical intersection dia- 
grams: 



r(i, 0, 0) = r^r^rj, = r^r^ = r^r 



(76) 



In order to compensate for identical products of subdia- 
grams, let us introduce the 



p : polynomial multiplicity. 



(77) 



For the current product of three vertex functions, the 
combinatorial factor p(n l7 n 2 , n 3 ) for r(ni, n 2 , n 3 ) derives 
from the binomial coefficient of the generating function: 



E 

k 1: k 2 ...=0 



3 

ki, k 2 , k 3 , 



reducing to the three representative cases of the index- 
vector k = (fel, k 2 , • • -) : 



(1=1,1) 
(2,1,0) 
(3,0,0) 



n 1 =n 2 = n 3 



p = 6 
P = 3 
p=l 



Taking this into account and the contraction multi- 
plicity m defined in (48), the generating function $1.3 



reduces to the product of three vertex functions: 

cr(r fc+3 ) 



$i,s(T) = 



fe=0 



(fc + 3) 



fc+3 



(79) 



E E 



a(n 1 ,n 2 ,n 3 ) 



711,12, ".3 = 1 »1,»2,'3 

m(ni, 722,^3), 1 ~ ni 



(ni + n 2 + n 3 + 3)! 

1 



ni \n 2 \n z \ (80) 
1 



p(rai, 712,713) ni! 



(— r ni v - (— r™ 2 v - (— — r" 3 v - 



r>2 



713 



I x C >W 



= *Y1 V(A) tlt2 V(B) t2t3 V(C) l3ll 



(81) 



with the constant k to be derived in the following. 

This result summarizes the central idea of the current 
work and requires some commends. First, the structure 
of ( 79 ) is uniquely defined by the virial expansion of the 



free energy |18|) and the requirement that its Euler form, 
suitably multiplied by powers of p and integrated over 
the particle coordinates, yields the free-energy functional 
density (63). In the next line (80 1, the representation 



( 64 1 of the triangular diagrams has been inserted, with 
only the indices of the backbone diagram ii,i%,i 3 ex- 



plicitely written out. After exchanging the sums in (81 1, 



the vertex functions are a convenient abbreviation for the 
resummed subdiagrams 



V ili2 (A) 



E 



1 



(82) 



The definition of the vertex function agrees with ( |61[ ) 
and introduces the factorials in ( 80 1 necessary to replace 
cr by <7. The prefactor k = am/p is therefore a function 
011 the equivalence classes of triangular diagrams and can 
be evaluated by inserting (68), (78), and (49). The de- 
tailed calculation can be found in appendixjX and yields 
the overall constant: 

K(n 1 ,n 2 ,n 3 ) = - for all r(ni,n 2 ,n 3 ) , (83) 
6 

independent of the particle numbers or equivalence 
classes. Instead, the prefactor corresponds to its coef- 
ficient of the backbone diagram k(0, 0, 0). 

The splitting of intersection diagrams into subdia- 
grams restricts the possible dependence of k on the num- 
ber of external particle-lines. Because of the invariance 
of the intersection centers under permutations, the pref- 
actor is either a function of the sum of particle num- 
bers /(ni +n 2 + 713), its product g{ni)g{n 2 )g{n 3 ) or the 
product of both. For a general diagram, whose backbone 
diagram of pairwise intersecting particles has k intersec- 
tion centers T(ni, . . . ,rik), this generalizes to functions 
invariant under the automorphism group 



( 78 ) K(m,...,n k ): gK = K for g € Aut(f (0, . . . , 0)) , (84) 



with r(0, . . . , 0) as the intersection diagram with all ex- 
ternal particle-lines removed. This indicates, but not 
proves, that the free-energy contribution for a given back- 
bone diagram can be represented by vertex functions 
that only depend on the number of paired particle-lines. 
The prefactor for all diagrams coincide and correspond 
to f(0,...,0). 
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v v h 



Vi 



11. ..Ik 



A A A. 



«1 %2 «1 »2 IS «1 «2 »* 

FIG. 9. Graphical illustration of fc-vertex functions: The 
crossed circle represents the resummation of external, un- 
paired particle-lines located at a single intersection center. 
Whereas the outgoing lines correspond to internal particle- 
lines of the backbone diagram, each carrying one of the fc 
particle indices. 



Translating back the diagrammatic representation ( 82 ) 



of the vertex function to the density dependent func- 
tional, we define the "vertex function" of k internal par- 
ticle lines 



oo 



(85) 



/ * — ' m + n 

J m=0 



lii...i n ki...k m )\ j 

)Pkx ■■■Pk m dn v 



as the sum of Eulcr forms over prime subdiagrams. Here 
we also used the observation of ( 62 ) that the integration 
over n v factorizes for each intersection center. Each k- 
vertex corresponds to a resummed intersection center, 
indicated by a crossed circle in Fig. [9j and k internal 
particle- lines. 

In terms of vertex functions, the Roscnfcld functional 
is written as the 0-vertex 



(86) 



= / V(r a ) d A r a 



whereas the Euler form of ( 79 1 translates into the product 
of 2- vert ices 



4>iAr a ,r b ,r c ) := - J V™ {? a )V l ™ {n)V^ (»*) 

Vertex functions are the building blocks of FMT func- 
tionals and, in analogy to the derivation of the 0-loop 
order ( 25 ) , can be written as a function of n v . Setting 



k = 0, Eq. (851 reproduces the Rosenfeld functional: 



J k=l 

/oo 
"xK)* -1 + (* - 1) C a ^n ai n a2 n k v - 
1 1 



+ (k — l)(Jfc - 2) C aia2a3 n ai n a2 n a3 n^ 3 dn v 

«1«2, 



n v 



l-n v 



2C a 



1 



(1 - n v y 

which yields the 0-vertex function: 

V(r a ) = -n x In (1 - n v ) + C aia2 n ai n c 



dn v , 



1 - n v 



''Ql ' V QL2 Ct'3 



1 

(l-n„)2 



(89) 



Correspondingly, the 1-vertex is a function with a sin- 
gle internal particle-line. However, as the virial expan- 
sion of the free-energy of hard particles only depends on 
star-graphs, no such term will occur. Nonetheless, the 
1-vertex provides the first correction in the Mayer ex- 
pansion of soft potentials with a hard-body center. We 
will therefore note its form for completeness: 



V^(r a ) = -w£ ln(l - n v ) + C a ^u l a \n a2 - 

1 



1 



^"10203, ,«i „ „ 

w ax "a2 "03 (-{ „ \2 



(1 - n v y 



(90) 



The analogous calculation with two fixed particle-lines 
and pi 1 , pi 2 removed 

/oo 1 
K{Yf-^)p i3 ...p ih dn v 

k=2 

/oo 
$>-l)C A - 
t — o 



iA 2 , Ai , .12 „k-2 



k=2 



+ (k - l)(fc - 2) C QlQ2tt3 «n Q3 ^- 3 dn v (91) 



_ / ^A 1 A 2 , Ai , Ai 
- ' ° u; A 1 UJ A 



1 (1 - 



w Ql w a2 n Q3 



yields the 2-vertex function 



t/«iZ2/V "\ — riA\A2, Ax . Ai 
V 7 (Taj - O W Ai W A2 j 



(1-n,) 3 



1 



(92) 



(l-n„) 2 



For vertices with more than two particle-lines, the sum 
reduces to one term only 

/oo -, ^ 

/OO 
7 



k—n 



v , .('1, Ai , A3 , ,*4 , ,in)„fe— n J„ 
X lv A 1 u 'A2 UJ A 3 uj v ■•■W v n v dn v 



(93) 



{n-l)\C AlMM 



v , ,»2 , ,*3 , ,14 , An 



dn v 



(1 - 

whose integration yields the vertex function 

V h - in (r a ) = (n-2)\C AlA2A3 (94) 



V , ,*2 , A3 , Ai , An 

X UJ A 1 u} A 2 U} A 3 UJ v ■ ■ ■ U v 



(1-n,)™- 1 



for n > 3 
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symmetrizised in the particle indices i\, . . . , i n . 

An important result of the resummation process is the 
pole structure of the generating functions. The vertex 
with k > 2 particle lines has a pole at least of order 
k — 1 at packing fraction n v = 1. This shows that the 
influence of diagrams rapidly decreases with their number 
of intersection centers. For the free-energy functional of 
hard spheres, this explains the success of the Rosenfeld 
functional. The leading correction <f>\^ will then be of 
order —3 and only take affect at high densities or strong 
angular correlations between particles, i.e. the solid state 
of the statistical system. 



V. THE SUBTRACTION SCHEME AND FIRST 
ORDER CORRECTION 

Using the representation of vertex functions simplifies 
the derivation of new functionals <S> g ,h for a given back- 
bone diagram. However, as each order of the expansion 
approximates an infinite subset of Mayer diagrams, it 
is not possible to simply add its contributions. For ex- 
ample, the Rosenfeld functional $q,i approximates the 
third virial contribution by its contracted form, whereas 
$1.3 contains its exact integral. The naive sum of both 
term $ ,i + $1,3 therefore causes a double counting of 
diagrams. How to compensate such terms by subtrac- 
tion will be shown in |V A[ followed by the resummation 
of Mayer ring diagrams and subsequent comparison with 
the White Bear II functional in IV Bl 



A. The first order correction to Rosenfeld's 
functional 

The splitting of intersection diagrams into subdia- 
grams and their subsequent resummation is a local map- 
ping, reflected in vertex functions that each depend on 
one intersection center only. Thus the global information 
about the topology of the original Mayer diagram is par- 
tially lost. Taking Fig. [7] as an example, the root points 
of different subdiagrams, indicated by i,j,k, are statisti- 
cally independent, and their Euler form has thus to van- 
ish for concurrent intersection centers. This, however, is 



in contrast to the functional (87 1, which is none-zero in 
the limit of coincident intersection centers B,C —> A. 

In order to obtain a physically consistent result, the de- 
generate contribution for r a = f% = r c has to be removed 
from the integral: 



1,3 



<pi,3(ra,rb,r c ) 
x d'yadjbd'jc 



B,C— >A 



(95) 



The subtraction of the degenerate part from the func- 
tional solves two problems: First, it removes the con- 
tractions incompatible with the Mayer diagrams. But, 
as a second effect, it also removes the free-energy contri- 
butions of lower diagrammatic orders, that correspond to 



identical Mayer clusters, but at different order of approx- 
imation. For example, the leading term of <^i,3(r a , f%, f c ) 
is the exact third virial integral, whereas the 3-particle 
contribution of the Rosenfeld functional Fq_i contains 
only its approximated form as a contracted diagram. 
Subtracting the contracted diagrams in ( 95 1 removes 



therefore its corresponding term in the 0-loop order of 
Fq,i + Fi,3- 

Removing unphysical terms from the free-energy is a 
common step in the regularization of loop integrals in 
quantum field theory. Because of this formal resem- 
blance, we will call the current subtraction scheme the 
"regularization" of intersection diagrams. 

Using the intersection algebra and the contraction 
rules for weight functions, the consistency of ( 95 1 under 



regularization can be shown explicitly. Leaving out the 
particle densities in (87), the contraction of the leading 



contribution of the 2-vertex product 



1 {c a ^lj Ax u a% y _ 1 



J x 



6 (l-n v ) 3 6" 

3-2 V 



reproduces the contracted 3-particlc virial and cancels its 
corresponding part in ( 25 ) . In summary, the functional 



Fr — Fi 



0,1 



-Pi, 3 



(97) 



is exact up to the third virial order and approximates 
all triangular and completely connected Mayer diagrams. 
The explicit dependence on 2-point densities now resolves 
the artificial degeneracy of Fq 1 in the orientational de- 
grees of freedom. Furthermore, the triangular diagrams 
describe distance correlations beyond the hard-particle 
diameter and therefore exceed the Percus-Yevick approx- 
imation. 



B. Approximations and the White Bear functional 

With the given set of rules, the derivation of functional 
corrections of arbitrary order becomes possible. But, as 
is often the case, including higher order terms does not 
guarantee higher orders of precision. Possible reasons are 
the divergence of series and the increasing calculational 
efforts to evaluate higher order terms. Actually, the ben- 
efit of the Rosenfeld functional and the vertex functions 
is their growing order in 1/(1 — n v ), which ensures a fast 
converging expansion away from its singularity. More re- 
strictive is therefore the second aspect, how to evaluate 
and minimize higher order terms, whose three and more 
intersection centers define non-local functionals. 

The mathematical framework necessary to evaluate 
such ring diagrams has been developed by Wertheim and 
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FIG. 10. $1 is the generating functional for all Hoop Mayer 
diagrams, whose contraction provides a first order approxi- 
mation in n v to the Rosenfeld functional. 



applied to the third virial integral for ellipsoidal geome- 
tries at constant particle density [32 35j . In the explored 
range of aspect ration L < 10, Wertheim found excellent 
agreement with results of computer simulations. 

The same mathematical methods apply to the func- 
tional (97), with the tensorial products of the normal 



vectors in ( 10 1 represented by spherical harmonic func- 



tions and the convolute of weight functions decoupled 
by a Radon transformation. The latter reduces to a 
Fourier transformation in the case of coinciding intersec- 
tion centers. The subtracted part of (95 ) can therefore be 



evaluated in the same way as the Rosenfeld functional. 
Nonetheless, evaluating and minimizing the functional is 
still a complicated mathematical problem, wherefore the 
development of efficient approximation strategies will be 
an important future goal. 

One possible ansatz is the improvement of the analyt- 
ical form of the n^-dependence of the functional. The 
previous results for the 0-loop order and the vertex func- 
tions suggest an expansion in powers of 1/(1 — n v ). This, 
however, is a result of the chosen resummation strategy. 
It is important to observe that different selections of di- 
agrams will also yield a different analytical structure in 
n v . As has been discussed in |25| , what marks the vertex 
functions as special is their possibility to combine with 
any alternative resummation scheme because of the fac- 
torization of Mayer diagrams into completely connected 
subdiagrams. 

In order to improve the n^-dependence of the 0-loop 
order, it is necessary to go beyond the starfish graphs. 
The first choice is therefore the set of 1-loop diagrams 



h=3 



(98) 



As a further approximation we restrict the functionals 
^g,h to their backbone diagrams, which results in the 



scries shown in the first line of Fig. 10 



The functional $i for Mayer ring-diagrams has already 
been derived in [25] as the generating function of Mayer 
bonds fi r i 2 , written in the matrix notation 



M lll2 := C A ^^yi Pil 



(99) 



Together with the symmetry factor cr(fc) = I /(2k) for a 
ring of k particles, the sum can be rewritten in closed 



form 



OO 



k=3 



%2i3 ■ ■ ■ Mi k i t 



(100) 



■- In (1 - M) - -M nn - -M ili2 M l2ll 



with a logarithmic singularity at M = 1. It is tempting 
to assume that this divergency corrects the pole n v — 1 
of highest packing fraction of the 0-loop functional to a 
physically realistic value that depends on the geometry 
of the particles. However, minimizing $i will be even 
more ambitious than that of Fk> 

As we are only interested in the n„-corrections of the 
Rosenfeld functional, it is sufficient to contract (100) to 
one intersection center. Using the notation 



r, 



— < F F 



(101) 



to indicate the contraction of a Mayer diagram, the log- 
arithmic part of (100) can be expanded in orders of n v : 



I 1 

--<ln(l-M) >=£ Tk <M h > 

k=l 

°° l 1 

E u K{ k Tii - ik)pk dnv 
fe=i 

/oo 1 „1 



d(tn v ) (102) 



fc=i 

1 [ "x 

2 / 1 - tn 



2C 



^Ql ^Q'q ^tt3 



d(tn v ) , 



where the integration over the scaling parameter p — > tp 
absorbs one factor of 1/fc. Evaluating the integral yields 
the logarithmic contribution 



- - < ln(l-M) >= --n x ln(l-n„) 



1 



4n„(l - n v ) 



(103) 

[2(1 - rc„)ln(l - n v ) +2n v ] 



2n2(l - nv y 



- [2(1 - n,,) 2 In (1 - n v ) + 2n v - 3n 2 v 



The remaining two terms of ( 100 ) are the traces of M 
and M 2 , which are functions of one intersection center 
only. Thus, using the subtraction scheme introduced in 
the last paragraph, both contributions would vanish af- 
ter regularizing the 1-loop diagrams. Unfortunately, this 
mechanism does not apply for the current approximation 
and a better understanding of these two terms is neces- 
sary. 

Deriving the trace of the contracted form of M 



< M hil > =< C 



MA 2 ,M ,M 



> Pil 



(104) 
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yields the weight density of the Euler characteristic. It 
therefore removes in ( 100 1 the case of only one intersect- 



ing particle. This geometric interpretation is consistent 
with Fig. [lOj where at least three particles have to in- 
teract at one intersection point. Generalizing this obser- 
vation to the calculation (102 1, we have to subtract the 



contributions of K(H n D k ~ L ) from its sum. The deriva- 
tion is analogous to the calculation of the 1-vertex (901, 
only with the symmetry factor replaced by cr(l) = 1/2 
and the sum restricted to Ed D k ~ 1 : 



If 1 

- o / E C XV - V n x n k v ~ 1 dn v = -n x In (1 - n v ) (105) 
fc=i 

The calculation for the second term M 2 is similar but 
allows two possible contractions: 

< M ili2 M l2ll > 

— < O UJ Ai UJ A2 U UJ A 2 UJ A 1 > PhPi2 

< M lll3 M l3l2 > Pl3 (106) 

— < O W Ai W A3 0J A 3 UJ A 2 > PtiPiiPi-i 
_ t ^a 1 a 3 a 2 

corresponding to intersections of E 2 and E 3 . But as both 
contributions derive from the same order in M, they nec- 
essarily have to follow from the same C-matrix. With the 
symm etry factor <j(2) = 1/4, the terms to be removed 
from ((T03| are the Euler forms of K(Y, 2 n D^ 1 ) and 



if (E 3 fl D k x ) multiplied by the same order of D k as 



in ( 105 ). The calculation therefore parallels ( 92 ) with n 



removed: 



1 1 7? 

^ V PiiPi 2 n v — n ai n a . 2 - — — 

— o n- ai n a2 n a3 j- r» ■ 
4 (l-raj 2 



(107) 



The new functional is the sum of the three contribu- 



tions (103), (105), (107) and the 0-loop order (25) 



$ = -n x ln(l -n v )(l + (t) (1) ) 
+ C a ^n ai n a2 —^—{\ + ^) 



1 - n v 



(108) 



£J a l a 



2«'3 



1 



(1 - n v y 

with the three correction terms: 



(l + (3) ) 



= 



1 



4 2) = 7^( 2 (! - nv) In (1 - n„) + 2n„ - rt 2 ) 



(109) 



v K - ( 2 ( X _ Tl «) ln(l - n^) +2n„ - 3n„ + • 

In this form, the result can be compared to the White 
Bear II functional, introduced in [2TJ [55] , which combines 



the Boublik-Mansoori-Carnahan-Starling-Leland equa- 
tion of state [40] with the structure of the free energy 
functional, determined by the scaled-particle differential 
equation [3], It is therefore not a purely geometrically 
motivated approach as Rosenfeld's FMT, but takes into 
account the numerically derived virial coefficients up to 
the eighth's order, combined in a generating function |41) . 
The WBII functional is independent of positional cor- 



relations and therefore of the same structure as ( 108 ) 
with the corresponding correction terms: 



6 (1) 







*WB = ^[2(l-n„)ln(l-n„) + 2n l ,-n 2 ] I I 10) 



4 WB = ^[2(l-n„) 2 ln(l-n„) + 2n„-3n 2 +2n 3 ]) . 



The terms </>k, </>wb are similar, deviating only in their 
numerical prefactors and the n 3 -term of </>( 3 '. Compar- 

we find excellent 



(2) 

ing the curve shapes of (fr K and 



6 (2) 



agreement despite this property. Whereas the different 
n 3 -terms of cause a significant change in the curve's 
curvature. Nonetheless, it is remarkable that the analyt- 



ical terms of the infinite sum of Mayer diagrams ( 103 ) 
are in good agreement with their analytical counterparts 
in the White Bear II functional. 



VI. DISCUSSION AND CONCLUSION 

The current article has shifted the previous perspective 
of [35] from differential geometry to the algebraic rules 
of the Euler form. It systematically generalizes the FMT 
functional from Rosenfeld's 0-loop order to any number 
of intersection centers and develops several new math- 
ematical tools for the efficient manipulation of weight 
functions. 

It has been shown that Mayer's star graphs have a 
uniquely defined dual representation in intersection di- 
agrams. These allow an intuitive picture of the Euler 
form and its decomposition into weight functions. Using 
this graphical description, we developed the contraction 
method as an approximation of the underlying Mayer di- 
agrams. 

Removing the external particle-lines from an intersec- 
tion diagram defines its backbone graph, which corre- 
sponds to a unique contribution $g t h in the expansion 
of the free energy functional. The successive attachment 
of external particle-lines to its intersection centers pro- 
vides a resummation process that significantly improves 
the original virial expansion in particle densities. The re- 
sulting vertex functions then replace the Mayer functions 
as the building blocks of the FMT functional. 

Resummation is an essential step in the derivation of 
the functional as it generates the pole structure 1/(1— n v ) 
of the packing fraction and yields a generic convergence 
criterion for the expansion in intersection centers. How- 
ever, its factorization in vertex functions could not be 
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proven in general. The problem relies in the symmetry 
factors of the infinite sum of Mayer diagrams. For the ex- 
emplary case of triangular graphs it could be shown that 
such a splitting indeed exists as the free-energy prefactor 
for all such diagrams agrees. 

In the current case of triangular diagrams, we first de- 
termined the invariance groups under labeling which de- 
fine four equivalence classes. Furthermore, taking into 
account the contraction and polynomial multiplicities of 
their individual vertex contributions, it could be shown 
that each triangular diagram contributes the same nu- 
merical prefactor. This allowed the simple factoriza- 
tion into fully contracted subdiagrams for any triangular 
graph. For general diagrams, however, we were only able 
to show that the prefactors for a given backbone dia- 
gram transform under the same automorphism group as 
the diagram itself. Nevertheless, this is enough to de- 
duce that the resummed free energy functional for any 
Mayer ring-diagram likewise factorizes, with a common 
prefactor that only depends on its number of intersection 
centers. But for more general functional contributions of 
two and more loops, it would be an important step in 
our understanding to obtain a general expression for this 
structure. 

The resummation of diagrams also involved the regu- 
larization of the functional. This process removes those 
intersection terms that are either incompatible with the 
Mayer diagrams or contributions of lower order function- 
als which are now replaced by terms of lesser approxima- 
tion. Most terms of the subtraction scheme are of zero 
measure and thus can be ignored completely. A prac- 
tical application is the summation over all completely 
contracted Mayer ring-diagrams. The analytical form re- 
produces the White Bear II functional to good accuracy. 
However, it is the subtraction of the irregular parts of first 
and second order that cause the deviation in that does 
not match with the White Bear result. If this discrep- 
ancy could be clarified it would be possible to derive even 
higher order corrections and exceed the current precision 
of numerically obtained functionals. 

The White Bear II functional provides only a small 
correction to the bulk properties of the Rosenfeld func- 
tional. This suggests that the improvement of the an- 
alytical structure in n v will be less important than the 
inclusion of terms which resolve the orientational degen- 
eration of the 0-loop order and to go beyond the Percus- 
Yevick approximation. Both deficits are resolved by the 
three-center term derived in this article. It will there- 
fore be an important next step to find efficient numerical 
methods to minimize Fx- 

The new objects entering this functional are the 2- 
point densities. Whereas the 1-point functions depend 
on a single vector field moving over the surface of a par- 
ticle, the 2-point functions determine the correlation be- 
tween vector fields at two different particle points. Thus, 
1-point functions reproduce the "classical" curvature de- 
pend information of the particle's geometry, whereas k- 
point functions define a completely new mathematical 



class of geometric invariant quantities. Fundamental 
measure theory introduces therefore new mathematical 
tools that not only give access to important physical 
problems of many-particle systems but also might pro- 
vide new answers to mathematical questions that cannot 
be covered by single vector fields alone, fc-point functions 
interpolate between the manageable but approximate de- 
scription of geometry as a tangential space and non-local 
geometric properties as, e.g., the maximal packing den- 
sity of particles in an embedding space. 

The derivation of a FMT functional starting from an 
infinite class of Mayer graphs, then translated to inter- 
section diagrams, and finally resummed into vertex func- 
tions is an inefficient approach. Given the diagrams' sys- 
tematic and our first experience with ring graphs, we ex- 
pect the existence of a simpler formulation that closely 
resembles a field theory with the vertex functions as its 
variabels. This would completely replace the dependence 
on Mayer diagrams and their symmetry factors and pro- 
vide a better understanding of the nature of FMT. 
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APPENDIX A 

In the following we will determine the automorphism 
groups and their characteristic parameters ( 65 1 , 1731) , and 



(831 for the triangular Mayer diagrams T(ni, U2, ns) of 
the type shown in Fig. [7| To the best of the authors 
knowledge, no general classification or systematic con- 
struction of these groups is known. We will therefore 
first derive the groups of the first seven diagrams and 
then deduce their generalization to all further cases. 

To put the formulation on a more formal level, consider 
any star diagram T and define a representation A of T by 
its labeling in Mayer's / functions: 



A 



prod(/y) 



(Al) 



By definition, any product of / functions is only uniquely 
defined up to their ordering 



^ ■ \fij fjii fijfkl fklfij} 

TrA(r) = A(r) . 



(A2) 



As any labeling of the nodes of a Mayer diagram is 
admissible, there exists n\ possible representations for an 
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4*^>»6 2 l^> ll 2 J ^^ l 3 



FIG. 11. The first Mayer diagrams of up to 6 particles with 
triangular substructure. (0,0,0), (1,0,0), (1,1,0), (2,0,0), 
(1, 1, 1), (2, 1,0), and (3,0,0) with the inner triangle defined 
by the numbers 1 — > 2 — s> 3. 



n-particle diagram, generated by operating with the sym- 
metric group S n on any representation A: 



£ S n : g\ = A' . 



(A3) 



For a discussion of the symmetric group see e.g. [12] , 
However, not all elements of S n generate a new repre- 
sentation. If the operation of g can be undone by a per- 
mutation 7T, it leaves the labeling invariant, denning a 
subgroup of S n - 

Aut(r) = {ge S n I g\(T) = tt\(T)} , (A4) 

so that 7T _1 o g = id. Clearly, the identity e € 5„ is 
element of Aut(r), and if 51,32 S Aut(r) so is <?r ,g% 
and any product of them. 

To simplify the notation, let us replace the / functions 
by square brackets: 



(A5) 



op- 



on which the permutation symbols (1, 2, 3, . . .) of S„ 
erate by cyclic permutation of particle indices. 

As an example, consider the labeling of the diagram 



r(l, 1, 0), as shown in Fig. 11 with its Z 2 -symmetry 2 -H> 
3, 4o5: 



M-Ti.i.o) — /12/13/14/15/23/24/35 

= [12] [13] [14] [15] [23] [24] [35] 



(A6) 



The axial symmetry of the diagram only allows the 
identity and one further element as the automorphism 
group 



(1)(23)(45)([12][13][14][15][23][24][35] 

= [13] [12] [15] [14] [32] [35] [24] 



(A7) 



so that the order of Aut(r(l, 1, 0)) is 2. 

The graphical representation of the diagrams, shown 
in Fig. |11[ suggests a relation to point-groups. This is of 
course only true for the simple graphs under considera- 
tion, but simplifies the construction of the automorphism 
groups considerably. For up to six particles, the star di- 
agrams and group elements are as follows: 



TABLE I. The triangular Mayer diagrams for n < 6 particles 
are listed with their automorphism groups and corresponding 
characteristic numbers: the symmetry factors a, a, the poly- 
nomial p and contraction multiplicities m, and their resulting 
prefactors k for the free energy functional. 



r 


Aut(r) 


a 


a 


P 


m 




(0,0,0) 


s 3 


1 


1/6 


1 


1 


1/6 


(1,0,0) 


Z 2 x S 2 


6 


1/4 


3 


2 


1/6 


(1,1,0) 


E'i x z 2 


60 


1/2 


3 


1 


1/6 


(2,0,0) 


£1 x Z 2 x S 2 


30 


1/2 


3 


1 


1/6 


(1,1,1) 


G?6 


120 


1/6 


1 


1 


1/6 


(2,1,0) 


E{ x Si 


360 


1 


6 


1 


1/6 


(3,0,0) 


E1XS2X S 3 


60 


1/2 


3 


1 


1/6 



(1)(2)(3),(1)(23),(2)(13),(3)(12),(123),(132) 
(1,0,0) :Z 2 xS 2 

(1)(2)(3)(4),(1)(4)(23),(2)(3)(14),(14)(23) 

(1.1.0) :EfxZ 2 
(1)(2)(3)(4)(5),(1)(23)(45) 

(2,0,0) : E 1 x Z 2 x S 2 
(1)(2)(3)(4)(5),(1)(2)(3)(45),(1)(4)(5)(23), 
(1)(23)(45) 

(1.1.1) :G 6 

(1)(2)(3)(4)(5)(6), (1)(5)(23)(46), (2) (6) (13) (45), 
(3) (4) (12) (56), (123)(456), (132)(654) 

(2,1,0) : E\ x S 2 (A13) 
(1)(2)(3)(4)(5)(6),(1)(2)(3)(4)(56) 

(3,0,0) : E l x S 2 x S 3 (A14) 
(1)(2)(3)(4)(5)(6), (1)(2)(3)(4)(56), (1) (2) (3) (5) (46), 
(1)(2)(3)(6)(45), (1)(2)(3)(456), (1) (2) (3) (465), 
(1)(23)(4)(5)(6), (1)(23)(4)(56), (1) (23) (5) (46), 
(1)(23)(6)(45), (1)(23)(456), (1)(23)(465) , 



(A9) 



(A10) 



(AH) 



(A12) 



(0,0,0) : S 3 



(A8) 



where Ei is the identity element. The diagrams and 
groups are also listed in Tab. [T] 

These groups have already been identified by Riddell 
and published by Uhlenbeck and Ford [39] ■ Based on 
Polya's counting theorem 03 Hi], they also developed a 
counting formula [45 \ Hi , which determines the number 
of independent labelings of all Mayer diagrams for a given 
number of nodes and labels. But a corresponding formula 
for individual diagrams is still unknown. 

For the current case of triangular Mayer diagrams 
r(ni, 7i2, na), the automorphism groups can be derived 
from the seven cases ( A8||A14 ). By comparing their 
groups with the diagrams of Fig. [TT] one observes that 
the three attached completely connected subdiagrams 
contribute the symmetric group S ni x S n2 x S n3 ■ For 
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TABLE II. Complete list of the triangular Mayer diagrams 
and their automorphism groups, their intersection symmetry 
numbers a, polynomial multiplicities p, and free energy coef- 
ficients K. 



V(n 1 ,n 2 ,n3) 


Aut(r) 


a 


P 


K 


m / n-2 n 3 


El X Sni X S n -2 ^3 


1 


6 


1/6 


ni = ni n 3 


Ei x Z 2 x 5^ x S„ 3 


1/2 


3 


1/6 


ni = n 2 = n 3 


G3„+3 


1/6 


1 


1/6 


Til = 7l2 = 0, 7l3 = 1 


Z 2 x S 2 


1/4 


3 


1/6 



m 7^ n 2 7^ 713, the resulting automorphism group there- 
fore is £f x S ni x 5„ 2 x S'ng. For m = n 2 7^ n 3 , the 
diagram has an additional axial symmetry Z2, exchang- 



ing the two subgroups S„ 



And for m = 712 = TI3, 



the symmetry of the backbone diagram extends to the 
dihedral group D§ , whose semi-direct product with is 
denoted as G 3n + 3 . 

As an example, how to generalize the above seven 
cases, consider T(l, 1, 1) and its extension to T(n,n,ri). 
In the corresponding diagram of Fig. [TTJ one replaces the 
particles 4,5,6 by a completely connected subdiagram 
of n particles. The dihedral symmetry of the backbone 
graph remains unchanged under this operation, so that 
the group elements (A12) can be adjusted by the formal 



replacement of the particle indices 4, 5, 6 — > A, B, C with 
A,B,CeS n : 

(n,n,n) : G 3 „ + 3 : ( Al 5) 
(1)(2)(3)(A)(B)(C), (1)(23)(B)(AC), (2)(13)(C)(AB), 
(3)(12)(A)(BC),(123)(ABC), (132) (CB A) . 

The resulting automorphism group is thus the semi-direct 
product D e x S%. 

Comparing these invariance groups to Tab. [TJ the only 
diagram that drops out of this classification is r(l,0,0). 
Instead of Z 2 x E\, as expected, its automorphism group 
is Z2 x Si. The reason for this larger group is an addi- 
tional Z 2 symmetry of its diagram, exchanging the two 
subtriangles 1 — 2 — 3 and 2 — 3 — 4, using the number- 
ing of Fig. [TT] As has already been observed in Section 
|III[ the contraction of the completely connected subdia- 
grams is therefore not unique and we have to count the 
contraction multiplicity m — 2 for (1,0,0) and m = 1 for 
all other triangular diagrams. Adding the exceptional 
case (1,0,0) to the previous list of triangular diagrams, 
all automorphism groups have been identified and are 
summarized in Tab. |JlJ 

The irreducible set of labeled diagrams can now be 
obtained by operating with the subgroup 



S„/Aut(r r , 



(A16) 



on one representative element A(T„). The number of dif- 
ferently labeled diagrams is therefore the order of this 
group, conventionally noted by the symmetry factor (65) 



of the virial integral (17). For the contracted dual dia- 
gram r, the corresponding symmetry factor has been de- 
fined in ( 73 1 , which for triangular diagrams simplifies to 



o"(r(m, n 2 ,n 3 )) 



S ni x S n2 x S n 



Aut(r(ni,n 2 ,n 3 )) 



(A17) 



As has been discussed in section 



TV] splitting the 



virial integral into vertex functions (79), (80), (81 ) is not 
uniquely defined but yields a multiple counting of virial 
diagrams r(ni, n 2 , 77.3) by permutation of its polynomial 
factors into V ( ni ) lll2 T {n 2 ) l2ls T (n^ 311 . For k identical 
values of the triplet m, n 2 , n%, the polynomial multiplic- 
ity is 



P 



3! 



(A18) 



The overall symmetry factor k for the decoupled virial 
integral is therefore the product 



P 



(A19) 



whose numerical values are listed in Tab. U and TAB. ITfl 
The central result of the current discussion is the free- 
energy pref actor, which for all triangular diagrams has 
the same numerical value 



(A20) 



as necessary for the decoupling and resummation of the 
intersection diagrams into vertex functions. 

For more complex classes of Mayer diagrams, the 
derivation of the automorphic groups is similar. However, 
apart from the exceptional cases, it is simpler to replace 
the Mayer diagrams by "weighted intersection graphs". 
They follow from the maximally contracted intersection 
diagrams by noting the number of unpaired particle lines 
at each vertex center and the subsequent removal of those 
lines from the diagram, as is shown in Fig. [7] For planar 
graphs, as the triangular example, the invariance group of 
the Mayer diagram derives from the semi-direct product 
of the symmetric groups of the vertices with the invari- 
ance group of the weighted intersection diagram. As the 
reduced diagram can be drawn in the plane, the latter is a 
discrete subgroup of 0(2, R). For more general cases the 
weighted intersection diagram can always be embedded 
into a Riemannian surface T g of minimal genus g |36j . 
The invariance group of the weighted intersection dia- 
grams is therefore either a discrete subgroup of 0(3, R) 
for the sphere g = or a subgroup of Sp(2g, n) for a 
diagram embedded into a torus of genus g. This depen- 
dence of the automorphism group on the topology of the 
graph might give a first explanation why Polya's count- 
ing theorem is applicable only for diagrams of individual 
classes. 
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